Relativistic and correlation effects in atoms 



B. P. Das, K. V. P. Latha, Bijaya K. Sahoo, Chiranjib Sur, Rajat K Chaudhuri 
Atomic and Molecular Physics Group, Indian Institute of Astrophysics, 

Bangalore - 560 034, India 
and 
D. Mukherjee 

Indian Association for the Cultivation of Science, Kolkata - 700 032, India 

2nd February 2008 

Abstract 

This review article deals with some case studies of relativistic and correlation effects in atomic 
systems. After a brief introduction to relativistic many-electron theory, a number of applications 
ranging from correlation energy to parity non-conservation in atoms are considered. There is a special 
emphasis on relativistic coupled-cluster theory as most of the results presented here are based on it. 
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1 Introduction 

One of the most important milestones in the development of theoretical atomic physics has been the 
formulation and application of theories that can simultaneously treat relativistic and correlation effects 
in atoms. Following the early work of Swirles on relativistic Hartree-Fock or Dirac-Fock (DF) theory 
Grant made pioneering contributions to the numerical and angular momentum aspects of this theory j2] 
which paved the way for further advances in the field. Multi-configuration Dirac-Fock (MCDF) 00] and 
relativistic many-body perturbation theory (RMBPT) codes were developed in the mid 1970s and 
early 1980s calculations based on them soon followed 1^1 El 13 H] . During the 1980s and 1990s, these two 
theories were appHed to a wide range of atoms and ions to study a variety of properties T0i, |nl El Ej ■ 
The extension of coupled-cluster theory to the relativistic regime during the last decade is indeed a very 
significant development [ll|. Linear and non- linear versions of this theory have been successfully used 
in performing high precision calculations of a number of different atomic properties ^15, 16, 17, l^j. 

The present review is by no means comprehensive; it mainly highlights some of the work on relativis- 
tic and correlation effects in atoms undertaken in our group. Unlike molecules, a number of different 
relativistic many-body calculations have been carried out on atoms using a variety of methods. Rela- 
tivistic many-body calculations on atoms are currently much more advanced than those on molecules 
|19[ I2()[ f?Tl 122) . In addition to the inclusion of the Breit interaction, certain types of QED effects have 
also been included in atomic calculations. It will take several years before molecular calculations reach 
this level of sophistication. Relativistic many-body calculations of parity and time reversal violations in 
some atoms have been performed to an accuracy of better than 1%. These calculations in combination 
with accurate experiments are now poised to test the Standard Model (SM) of particle physics. It is not 
clear at the present time whether it would be possible to achieve something comparable from studies on 
discrete symmetry violations of molecules. 

The organization of the paper is as follows : Section 12 deals with the Dirac-Coulomb approximation 
and the section following it (section O , touches upon the Breit interaction and QED effects. Section 
is an overview of relativistic coupled-cluster theory which has been used in the majority of the cal- 
culations considered here. In sectional we present the basic ideas underlying two physical effects that 
are relativistic in origin - fine-structure splitting and permanent electric dipole moment of atoms arising 
from the electric dipole moment of an electron and have given the results of some representative cal- 
culations. The enhancement of relativistic effects in heavy atoms along with the influence of electron 
correlation is discussed in sectional with reference to correlation energy, hyperfine interactions and parity 
non-conservation in atoms. In the last section we make some concluding remarks. 
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2 The Dirac-Coulomb Approximation 



For an A^-electron atom, the relativistic Hamiltonian is given by 

o2 



JV N 

H = J2[cai■P^ + f3mc^ + VN{ri)]+^ — , (1) 



where a and /? are given by cti = ^ ^ J and = ^ q j ^ Pauh matrices and / 

represents the unit matrix. VAr(ri) is the nuclear potential at the site of the zth electron and the last 
term is the Coulomb interaction between the electrons. H defined above is known as Dirac-Coulomb 
Hamiltonian which is clearly not covariant. 
This Hamiltonian can also be written as 

H = Y^ho{^)+J2f. (2) 

The electron-electron interaction can be approximated by an average potential where each electron 
moves independently in an average field caused by the nucleus and the other electron. This is the 
independent particle model which is the starting point of most atomic physics calculations. This can be 
put into a mathematical footing by partitioning the full Hamiltonian in the following way : 

H^Ho + Ves, (3) 

where 

N 
i 

is a sum of the one electron operators, 

hQ{i) = cai ■ Pi + (37710^ + U{ri). (5) 
It is customary to assume U as the Dirac-Fock potential and 

can be treated as a perturbation if there are no strongly interacting configurations in the system. The 
many-body atomic state |5'(r, J, M)) is an eigen function of the Dirac-Coulomb Hamiltonian and satisfies 
the equation, 

H\^ir,J,M))^E\^{T,J,M)), (7) 

where J, M are the total angular momentum quantum numbers and T is the quantum number which dis- 
tinguishes each of the atomic states. These states are expanded in terms of the determinantal wavefunc- 
tions which in turn are built from the single particle orbitals. If |$(r, J, M))'s denote the determinantal 
wavefunctions, then, 

\^iT,J.,M))=Y,Ck\M^,J,M)). (8) 

k 

The coefficients CfeS are determined by the choice of the theory. The single particle orbitals are the 
two-component Dirac spinors. 



1^"-) = I f .n"i1f " ), (9) 

Km 



r \ iQnK{r)x 



where n and m are the principal quantum number and magnetic quantum number respectively, k is a 
quantum number given by 

I fOTj = l-- 

-{l + l) forj-/ + i ' 
where I is the orbital angular momentum and j is the total angular momentum of an electron. 
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Figure 1: Positive and negative energy states 



The solutions of the Dirac equation admit both positive and negative energy states [22l and this is 
shown in figure For a free particle, only continuum states exist above mc^ and below —mc?. However, 
electrons in an atom that are acted on by a relativistic mean-field potential in addition to continuum 
states above mc^ and below — mc^bound states do exist in the interval —mc^ and mc^. The variational 
principle fails due to the presence of negative energy states ^23]. The radial parts of the large and small 
components of the Dirac spinor are expanded in terms of Gaussian functions [2^ as follows 

i 

QUr) ^Y^^LgUr). (11) 

i 

where S and L stands for small and large component respectively and the g's are Gaussian type functions 
of the form 

5^,M-C'i^'^""exp(-a.r2). (12) 

In the case of finite basis set expansions the condition of kinetic balance is appHed to prevent the 
variational collapse 26]. The kinetic balance condition [23|[23|2S| gives the relation between the large 
and small component of radial wave function as follows : 

5f.W = C^(^ + ^)5^.M. (13) 

In Eq. Ijl2|l and \l?t\ and Cff are the normalization factors for the large and small components 
respectively. 



3 Beyond the Dirac- Coulomb Approximation 

The electron-electron interaction can be treated relativistically, by including corrections to the Coulomb 
interaction. The leading relativistic correction to the Coulomb interaction is the Breit interaction [29] . 
where the interaction Hamiltonian is given by, 

^^^^^,^a^^ {a,-n,){a,-n,) ^ (14) 

Here the matrices ai,aj are built from the Dirac matrices and is the inter-electronic distance. The 
magnitude of the Breit interaction is smaller than that of the Coulomb interaction by a factor , where 
a is the fine structure constant and it can be included perturbatively or self consistently ^01. In addition 
to the Breit interaction, inclusion of QED effects like the self-energy and the vacuum polarization [HI] 
may be necessary for an accurate quantitative description of certain properties where relativistic effects 
axe important. 

The process involving the emission and absorption of a virtual photon by the same electron is known 
as self-energy. According to Dirac's theory, the vacuum consists of a homogeneous sea of negative-energy 
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Figure 2: Corrections to Coulomb interaction 

electrons. A bound electron in the atom can interact with an electron in the Dirac sea, thereby changing 
the charge distribution of the negative energy electrons compared to the free-field case. This results in 
the creation of electron-positron pairs and hence the vacuum behaves as a polarizable medium. This 
process is known as vacuum polarization. Figure|2lillustrates the self-energy and the vacuum polarization 
processes [32 . 

Only a few calculations of self energy and the vacuum polarization (which together give rise to the 
Lamb shift) have been performed on many electron atoms. As the Coulomb interaction due to the 
nucleus is much stronger than the electron-electron interactions in the inner shells of heavy atoms, it 
is reasonable to calculate the Lamb shift for such systems using the hydrogenic or screened hydrogenic 
approximation [22l- More sophisticated calculations of QED effects have been carried out in the past 
few years. The details of these calculations can be found in a review article by Shabaev |34j . 

4 Relativistic Coupled Cluster Theory : An overview 

In this section we will briefiy introduce the relativistic coupled cluster theory; one of the most pow- 
erful and accurate relativistic many-body theories. It is equivalant to all order relativistic many-body 
perturbation theory and has the virtue of being size-extensive . Most of the work described in the 
subsequent section are based on this theory. 

We start with the DF state |<i>) built out of four component orbitals given by Eq. Q, as the Fermi 
vacuum, and then the normal ordered Hamiltonian can be expressed as 

HN = H-{^\H\<i>)=H-EDF, (15) 

where H is the Dirac-Coulomb Hamiltonian. 

If we project ($|exp(— T) from the left we obtain the correlation energy (AE) and if we project 
any of the excited determinant ($*| exp(— T) we additionally get a set of equations which are used to 
obtain the T amplitudes. Using the normal ordered dressed Hamiltonian Hpf — exp(— r)iJAr exp(T) the 
corresponding equations for correlation energy and amplitudes become 

= AS, (16) 

and 

{<f*\HNm^O- (17) 

Here the state |<I>*) may be singly excited or double excited 1$^^) and so on. The indices a,b, - ■ ■ 
refer to holes and p, g, • • • to particles. We have considered the coupled cluster single and double (CCSD) 
approximation, where the cluster operator T is composed of one- and two-body excitation operators, i.e. 
T = T1+T2, and are expressed in second quantization form 

T = Ti + T2 = ^ [alaa] tl + \Y. H^'^baa) (18) 

ap abpq 

Contracting the ladder operators jSH] and rearranging the indices, the amplitude equations can be 
expressed in the form 

A + B(T)-T = 0, (19) 
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where A is a vector consisting of the matrix elements ($*| Hjsi |<i>) and T is the vector representing the 
excitation ampUtudes and BiT) is the matrix which depends on the cluster amplitudes so that Eq. Ijl9|l 

is solved self-consistently. For example, a typical contribution to the term HmT2 T2 is 

B^i - 1 E y^^rstz^i (20) 

dgrs 

Here Vdgrs is the two-electron Coulomb integral and is the cluster amplitude corresponding to a 
simultaneous excitation of two electrons from orbital a and d to p and r respectively. Diagrammatic 
techniques are used to obtain all the terms which contribute to this specific contribution. 

For an atom with one valance electron we first compute the correlations for the closed shell system, 
i.e. singly ionized atom using the closed shell coupled cluster approach. The reference state for the open 
shell system is 

K-'')^al\<P) (21) 

with the particle creation operator a^. Then by using the excitation operators for both the core and 
valance electron the exact state is defined as |37j : 

|M'f+i)=exp(r){exp(5,)}|$f+i). (22) 

Here {exp(S'fc)} is the normal ordered exponential representing the valance part of the wave operator. 
Here 

Sk = Sik + 5'2fe = E {'^l^k} Sfe + E {o'WgO.bak} sH , (23) 

k^p hpq 

where k stands for valance orbital. Sk contain the particle annihilation operator Uk, and because of 
the normal ordering it cannot be connected to any other valance electron excitation operator and so 
{exp(5fc)} automatically reduces to (1 + Sk). 
Then we can write the En. ll22|l as 

K+i) =exp(r)(l + 5fc)|$f+i), (24) 

and obtain a set of equations PE] 

(<i>f+l|i?^(l + 50|<i>r')=^e// (25) 



and 



Hn (1 + Sk) = Heff h7+' (1 + Sk) ■ (26) 



The Ea. (l26ll is non-linear in Sk because H^ff is itself a function of Sk. Hence, these equations have to 
solved self-consistently to determine the Sk amplitudes. 

The normahzed transition matrix element (i — > /) due to an operator O is given by 



A _ (>i'7+i|o|*f+i) 

_ (<t.7 + ^|{l + st}cxp(Tt)Ocxp(T){l+S.}|$f + ^) ^^'> 

^ ^(*7 + i|{l + St}exp(Tt)exp(T){l+Sj}|<I.7+i)« + i|{l+St}oxp(Tt)oxp(T){l + S.}K + i) ' 

whereas the expectation value of any operator O can be written as the normalized form with respect to 
the exact state |\I>^+i^ as 

^\ O 1*^^+1) _ {1 + 5t} exp(Tt)6exp(r) {1 + 5} 

^/ - ^7v+i) - {1 + 5t} exp(Tt) exp(r) {1 + S} ' 
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Table 1: Fine structure intervals for B-like ions in 



z 


Dirac-Fock (DF) 


Dirac-Coulomb (DC) 


Breit 


QED 


Total 


Experiment 


20 


37581.8 


38119 


-1616 


124 


36627 


36615(30) 


22 


57678.6 


58319 


-2241 


180 


56258 


56243(4) 


24 


84973.2 


85715 


-3015 


251 


82951 


829926(20) 


26 


121045.8 


121898 


-3958 


341 


118281 


118266(20) 


28 


167653.8 


168635 


-5092 


451 


163994 


163961(50) 


30 


226739.4 


227889 


-6437 


586 


222038 





5 Purely relativistic effects 
5.1 Fine-Structure splitting 

The fine-structure splitting is relativistic in origin, but is infiuenced by electron correlation. It occurs 
between the states having same values of the total orbital quantum number L, total spin quantum 
number S and different total angular momentum J. There have been many attempts to calculate this 
quantity for a variety of atoms in their ground and excited states jHH] ■ We present here the interesting 
case of the ground state fine structure splitting of boron which has been calculated by different relativistic 
approaches [39, 40. II. Boron is an open-shell atom with the configuration ls^2s^2p^. The relativistic 
configuration interaction method (CI) was used by Das et al jll]. The single particle orbitals used in 
the calculations there were obtained by the application of the variational principle. Consider an energy 
functional given by 

e ^yar{^r\H\^r) (29) 



where H is the Dirac-Coulomb Hamiltonian, |$r)is the rth configuration state function (CSF) and is 
given by 

2 Jr + 1 

Jr and Js being the total angular momenta of the rth and sth CSFs respectively. Minimization of e with 
respect to the single particle orbitals ^i, 

yields a set of differential equations which were solved self-consistently by using an appropriate numerical 
method ^2]. In this calculation, all relativistic configurations arising from ls^2s'^2p^, ls^2s^2p^and 
ls'^2p^ were considered. The Breit interaction and QED effects (self-energy and vacuum polarization in 
the hydrogenic approximation) were treated as first order perturbations. We give below the results of the 
calculation. The electron correlation contributions (difference of Dirac-Coulomb+Breit and Dirac-Fock) 
vary from -1078.8 to -5287.4 for Z = 20 to Z = 30. It is the evident from table[l]that the Breit 
interaction and the QED effects play an important role and their inclusion is critical in obtaining good 
agreement with experiments. 



5.2 Electric Dipole Moment of the electron 

The presence of a non-zero electric dipole moment (EDM) on a non-degenerate physical system would 
be a direct evidence of Parity (P) and Time-reversal (T) symmetry violations. An atom can have a non- 
zero EDM due to a non-zero EDM of it's constituent electron, under certain conditions. According to a 
theorem by Schiff, in 1963 |4^, the EDM of an atom vanishes even if it's constituents have non-vanishing 
EDMs. This theorem was based on the following assumptions : 

1. the constituents of the atoms are non-relativistic particles, 

2. the interactions between the particles in an atom are electrostatic, 

3. the EDM distribution of each atomic constituent coincides with its charge distribution. 
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By considering the relativistic effects in atoms, Sandars showed that an atom can have a non-zero EDM 
|44j . If an electron has a non-zero EDM de, the relativistic interaction of de with the internal electric 
field of the atom, is given by, 

i?/ = -^4/3^a^•^^^ (32) 

i 

where is the electric filed inside the atom; (li are Pauli matrices and /3 is Dirac matrix defined in 
section [2I This reduces to. 

Hi = -J2deff^■E^' (33) 

i 

in the non- relativistic Hmit. It is possible to express the relativistic form of Hj in terms of an effective 
one particle Hamiltonian, given by [3^ 

Hi = 2icPah} , (34) 

where c is the velocity of light, 75 = «7o7i7273 and 7i = I3ai. The Schrodinger equation for the 
unperturbed state ^) is 

i?o|*i°)) = 4°)|*i°^), (35) 

where \^^a') — exp(r(°))|$i°'') in coupled-cluster theory and i/o is the Dirac-Coulomb Hamiltonian. 

In the presence of EDM interaction, which is treated as a perturbation, the Schrodinger equation 
becomes 



H\m^) = , (36) 

where H ^ Hq + \Hi and l^-^) = exp(r(") + Ar(i))|$i°^). Here T('')and T^^^are the unperturbed and 
perturbed cluster amplitudes and the perturbation parameter, \ = de- The T^^^and T^^^ amplitudes are 
determined from the following equations : 



($*|i7jv|*) -0 (37) 



and 

The atomic EDM is given by. 



= 1 (38) 



for non-relativistic case. D is the electric dipole operator. The enhancement factor R is given by, 

""^^r ^WW) 

Following Coupled Cluster theory the equation (Eq. QHl reduces to 

($o|*o) 

An alternative Coupled Cluster approach to EDMs is given by Shukla et al ^^l- The values of de 
predicted by various models of particle physics are given in table 121 

The current best limit on the electron EDM comes from the Tl measurement (33 • The enhancement 
factor atomic thallium to the electron EDM R — —585, which is based on relativistic coupled-cluster 
calculation Comparing with experiment, the limit on the electron EDM is 

de < 1.6 X 10-2^e-cm. (42) 

The enhancement factor for atomic Cs (Z=55) has been obtained as i? = 130.5 from a method combining 
RMBPT and the MCDF approach gHI- The calculation done by Martensson et al fSD] gives R = 
114(1 ±0.03) for Cs. 
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Table 2: Value of de predicted by various models of particle physics 



Model 


de 


in e 


- cm 


Standard Model 




< 10 


-38 


Supersymmetric 


10- 


-26 _ 


10-28 


Multi-Higgs 


10' 


-26 _ 


10-28 


Left-right asymmetric 


10- 


-26 _ 


10-28 




Figure 3: Diagrams contributing to correlation energy 

6 Relativistic enhancements 
6.1 Correlation energy 

In the frame work of coupled-cluster theory, the expression for the correlation energy of an atom is given 
by, 

E,orr^mHN\^), (43) 

where Hn — e^^ Hj^e^ where H is the Dirac-Coulomb Hamiltonian described in section |21 

The diagrams given in figure contribute to the correlation energy, where the dotted Hues represent 
the coulomb interaction between the electrons, the solid Hue corresponds to the cluster operator T and 
the circle represents the one-electron operator. 

We have computed the above expression for the correlation energy using the coupled-cluster wave- 
functions for the closed-shell atoms, Xe^^, Yb'^^ and Hg^^ . The results are shown in tableEl 

From the above results it is clear that with the increase in the atomic number {Z), the relativistic 
effects become more prominent. The absolute magnitude of the Dirac-Fock contribution hence increases 
and that of the correlation energy decreases for xenon, ytterbium and mercury. 

The orbitals used in the calculation are expanded in terms of Gaussian functions of the type jfilj 

= r'=exp(-air2), (44) 

with fc = 0, 1, 2 • • • for s,p, d, • • • type functions, respectively. The exponents are determined by the even 
tempering condition |52) 



Table 3: Comparison of correlation and Dirac-Fock energy 



Atom 


Dirac Fock enerj 


;y 


^-^corr 




-0.74474960061E- 


h04 


-0.71694286411 




-0.14069217432E- 


h05 


-0.56956394691 


Hg80 


-0.19650686115E- 


h05 


-0.44792843655 



8 



Table 4: Details of the basis used in the calculation 



symmetry Total basis in No. of excited ag and /? used 

each symmetry orbitals 

7^2 13 8 0.00725 ; 2.725 

Pi/2 11 7 0.00755 ; 2.755 

P3/2 11 7 0.00755 ; 2.755 

^3/2 8 2 0.00775 ; 2.765 

^5/2 8 2 0.00775 ; 2.765 

/5/2 5 5 0.00780 ; 2.805 

/7/2 5 5 0.00780 ; 2.805 

97/2 3 3 0.00785 ; 2.825 

09/2 3 3 0.00785 ; 2.825 



a, =ao/9'"'- (45) 
The values of aoand (3 for different symmetries are given in tabled) 

6.2 Hyperfine interaction 

A nucleus may possess electromagnetic multipole moments, which can interact with the electromagnetic 
field produced by the electrons at the site of the nucleus . The interaction between various moments 
of the nucleus and the electrons of an atom are collectively known as hyperfine interactions [23 . This 
interaction produce shifts of the electronic energy levels which are usually much smaller than those 
corresponding to the fine structure splittings. 

The non-vanishing moments are the magnetic multipole moments for odd k and electric multipole 
moments for even k. The most important of these moments is the magnetic dipole moment (fc = 1) 
which is associated with the nuclear spin. The interaction of this particular moment with the electron is 
known as magnetic dipole hyperfine interaction. 

In general the hyperfine interaction is given by jj^l 

i7„/, = ^Af«.TW, (46) 

k 

where M^^) and T^^' are spherical tensors of rank k, which corresponds to nuclear and electronic parts 

of the interaction respectively. 

For the magnetic dipole hyperfine interaction 

where 7? is the Dirac matrix and Y^^ is the vector spherical harmonics. In Ea. H47ll the index j refers 
to the j-th electron of the atom and e is the magnitude of the electronic charge. The magnetic dipole 
hyperfine constant A is defined as 

^.,„fH^ (i||r<.)||J) 



I I ^J(J+1)(2J+1) 

where [ij^ is the nuclear Bohr magneton, [ij is the nuclear magnetic moment, J is the nuclear spin, J is 
the total electronic angular momentum. 

In En. lHTjl i'^^ is the single particle reduced matrix element of T^^^. The reductions of the single 
particle matrix element into angular factors and radial integral can be obtained by using the Wigner 
Eckart theorem. This single particle reduced matrix element is given by 

(«|| id) \\^) = - («|| {k + ^) f dr^InQj^L + QjlM , (49) 
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Figure 4: Goldstone diagrams for pair correlation (a,b) and core-polarization effects (c,d). Here a denotes 
a hole whereas v denotes valance orbital and p, 9, r... denote virtual orbitals (particles). The superscripts 
refer to the order of perturbation and the dashed Hues correspond to the Coulomb interaction. Particles 
and holes (labeled by a) are denoted by the lines directed upward and downward respectively. The 
double Hne represents the O (the hyperfine interaction operator) vertices. The valance (labeled by v) 
and virtual orbitals (labeled by p,q,r..) are depicted by double arrow and single arrow respectively, 
whereas the orbitals denoted by ® can either be valance or virtual. 



where (k|| C*^'^-' is the reduced matrix element of the Racah tensor and is equal to 

(-i)-'+^/v(2j+i)(2j'+i)(^ i ^,y{i,k,i'), 

with 

1 iil + k + l' even 
otherwise 

Here the single particle orbitals are expressed in terms of the Dirac spinors with Pi and Qi as large and 
small components respectively. 

In the calculation for Sa+we have used hybrid basis functions which are partly numerical and partly 
analytical . The analytical orbitals have the form of Eq. ijiiji . In table El the values of the magnetic 
dipole hyperfine constant (A) is given in MHz for ^^M^+LSgi and ^^'^Ba'^ for ground and one low 
lying excited state. In tableElwe have presented the contributions from Dirac- Fock (DF), pair correlation 
(PC) and core polarization (PC) effects. 

It can be seen from this table (table that for Mg^ the CP contribution is larger than the PC in 
magnitude for both the states. It is important to note that the former contribution includes the hyperfine 
interaction of all the core orbitals while only a specific valence orbital is involved in this interaction for the 



■K{l,k,l') 
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Table 5: Values of magnetic dipole hyperfine constant (A) in MHz for ^^Mg'^and ^^'^Ba'^ 



Atom Theory Experiment Others 

25Mc/+ States 

3si/2 592.86 596.25 602(8) gHI 

597.45 m 

3pi/2 101.70 103.4 (59] 

6si/2 4072.83 4018 4203.200 ^01 

6pi/2 736.98 742.04 



Table 6: Contribution of pair correlation (PC) and core polarization (CP) effect in magnetic dipole 
hyperfine constant (A) in MHz 

Atom Dirac-Fock PC CP 



25Mg+jSSI States 

3si/2 468.819 39.713 77.767 



137 



3pi/2 77.975 7.293 15.153 



6si/2 2929.41 663.20 465.91 
6pi/2 492.74 126.53 98.98 



latter (see Fig. However, the hyperfine constant A for i?a+ exhibits exactly the opposite behaviour. 
Even though Ba~^ has more core electrons than Mg^ , the relativistic enhancement of the valence (6s) 
magnetic dipole hyperfine interaction results in the value of PC exceeding that of CP. 

In tableHwe present the values of A for 6p^Pi/2 and 6p^P3/2 states of The Dirac-Fock values 

for these two states deviate from their experimental values in opposite direction, suggesting that the signs 
of the correlation contributions are opposite for the two cases. This is evident from the result of our 
second order relativistic many-body perturbation theory (RMBPT(2)) calculation. Electron correlation 
is dramatic in the case of the 6p^P3/2 state because of the large and negative core polarization (-840.6 
MHz). However, the value of A at this level differs from experiment by 48%. After carrying out a 
RCCSD(T) calculation this discrepancy reduces to less than 7% . The agreement of the ground state 
value of A with experiment is about 0.7%. These calculations highlight the power of the relativistic 
coupled-cluster theory to account for the interplay of relativistic and correlation effects in systems with 
strongly interacting configurations [62] . 

6.3 Parity non-conservation in atoms to neutral weak interaction 

The parity transformation can be expressed as — > — and the action of the parity operator P is 
given by 

P^iT) = i^i-T) , (50) 

where ipiT") is the wavefunction of a physical system. Parity conservation means that the system is 
invariant under parity transformation and the Hamiltonian H commutes with the parity operator, i.e. if 
Hp is the parity transformed Hamiltonian 
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Table 7: Magnetic dipole hyperfine constant for 6p states of P6+: a strongly interacting system. 
RMBPT(2) stands for second order RMBPT. Both RMBPT and RCCSD calculations are performed 
by our group. 



States 






Dirac-Fock 


11513.5 


918.4 


RMBPT(2) 


15722.5 


302.9 


RCCSD(T) 


12903.7 


623.2 


Experiment |^ 


13000 


583(21) 



Hp = PHP-^ = H 



and therefore 



H,P =0. (51) 



Hence it clearly implies that parity non-conservation (parity violation) means that its Hamiltonian does 
not commute with the parity operator P. 

Parity non-conservation (PNC) was discovered in the beta decay of ^°Co by Wu and co-workers in 
1957 following the prediction by Lee and Yang a year earher [OBI- This lack of mirror symmetry has now 
been observed in several systems and even in atoms which is an important phenomenon to study. The 
latest measurement of parity non-conservation in cesium with unprecedented accuracy (0.35%) has led 
to the discovery of the nuclear anapole moment |(i4j . 

The dominant contribution to PNC in atoms comes from the neutral weak current (NWC) interaction 
between the electron and the nucleus [65j- The effective Hamiltonian describing the interaction consists 
of two parts, one of which is nuclear spin independent (NSI) [SHI and the other is nuclear spin dependent 
(NSD) jSSIinZI- In this review article we will concentrate on NSI parity non-conservation in atoms. The 
NSI effective Hamiltonian is expressed as 

Hpr,c = ^QwY^ltP^^^)^ (52) 

with 



Qw = 2 [ZCip + NCm] ■ (53) 

Here Z and N are the number of protons and neutrons respectively and Cip and Ci„ are the vector 
(nucleon) - axial vector (electron) coupling coefficients whereas Gf is the Fermi coupling constant and 
p{re) is the normalized nucleon number density. The matrix element of Hpnc scales as and it has 
been treated as a perturbation. It is primarily because of this reason that the heavy atoms are considered 
to be the best candidates for PNC experiments. The total Hamiltonian is now represented by 

H = Ho + HpNc- (54) 

This perturbation causes the wavefunction to take the form \^) = + |^'^^), where and 

1^^^^) are the unperturbed and the perturbed part of the wave function respectively. 

The quantity that is measured in such an experiment depends on the interference of a parity non- 
conserving electric dipole transition amplitude {EIpnc) and an allowed transition amplitude corre- 
sponding to two atomic states of the same parity [0H|- From the theoretical point of view an accurate 
calculation of Elppfc must be based on a suitable and accurate relativistic many-body theory. In a 
recent review, Ginges and Flambaum have presented the current status of atomic PNC calculations 
and experiments. A number of many-body theories have been applied to calculate Elp^c matrix ele- 
ments. The results of these calculations in combination with the most accurate PNC experiment on Cs 
is in agreement with the Standard Model (SM) of particle physics 169) . 
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Table 8: ElpNc matrix elements for Cs and Ba 



Atom 


Transition 


ElpNci^eaoiQw/ - N)) 


Cs 




0.910 X 10"^^ 


Ba+ 




2.05 X 10"" 



We have formulated a new approach to PNC in atoms based on relativistic CC theory in an attempt to 
go beyond the existing calculations. In this formulation the excitation operators (both T and S) contain 
an unperturbed (superscript 0) and a perturbed part (superscript 1). For a single valence systems like 
cesium the wavefunction can be written as 



exp(T" + T(i)) {l + + ^(1)} |$fc) . 



(55) 



This equation follows from En. lf2ijl and can be derived easily [201 ■ The equations for determining T^^^and 
iS^^^ amplitudes are the following : 



l*o) 



PNC 



l$o) =0, 



(56) 



and 



{^t\HNSi'^ - A^o)^!) |$„) + i/^ {r(i) + Td)^")} + HpNc {l + Si"^} |$o) = . (57) 



The parity non-conserving electric dipole transition amplitude between atomic states 



) and 1^-/) = 



(0) 



^'y-^'' ) is given by 



{^f\D\M>, 



ElpNc = (5g) 

The preliminary results we have obtained using this approach are given in tabled These calculations 
have been carried out in the Dirac-Coulomb approximation with an universal Gaussian basis consisting 
of 13s, 12p, lid and 7/ function for Cs and 13p, 12p, lid and 8/ functions for Ba'^ . Our results for 
Cs is in reasonable agreement with Hnear relativistic CCSD(T) calculation |23 which yields EIpmc = 
0.0908(9) X 10^^^ in the same units as our calculation. 



7 Conclusion 

It is clear that considerable progress has been made during the past three decades on the relativistic 
many-body theory of atoms. However there are some open problems in this field. Perhaps the two areas 
that deserve most attention is the future are : 

(i) Relativistic multi-reference theories to treat a wide variety of open shell heavy atoms including 
rare-earths. Work in this area is in its infancy [72| . 

(ii) Incorporation of QED effects in a systematic way in the framework of relativistic many-body 
theory. 

One can indeed look forward to exciting new developments in relativistic electronic structure of atoms 
in the coming decade. 
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